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Abstract 

The lifetime splitting between the and -B^ mesons has recently been cal- 
culated in the next-to-leading order of QCD. These corrections are necessary for 
a reliable theoretical prediction, in particular for the meaningful use of hadronic 
matrix elements computed with lattice QCD. Using results from quenched lattice 
QCD we find r{B+)/T{B^) = 1.053 ± 0.016 ± 0.017, where the uncertainties 
from unquenching and 1 /nit corrections are not included. The lifetime difference 
of heavy baryons Hj^ and is also discussed. 
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Figure 1 : Weak annihilation (WA) and Pauli interference (PI) diagrams in the lead- 
ing order of QCD. They contribute to r(i3°) and r(i?+), respectively. The crosses 
represent |Ai3| = 1 operators, which are generated by the exchange of W bosons. 
CKM-suppressed contributions are not shown. 



1 Introduction 

In my talk I present work done in collaboration with Martin Beneke, Gerhard Buchalla, 
Christoph Greub and Alexander Lenz [|l]] . 

Twenty years ago the hosts of this conference showed that inclusive decay rates of 
hadrons containg a heavy quark can be computed from first principles of QCD. The 
Heavy Quark Expansion (HQE) technique [|[] exploits the heaviness of the bottom (or 
charm) quark compared to the fundamental QCD scale Aqcd- In order to study the 
lifetime of some 6-flavored hadron H containing a single heavy quark one needs to 
compute its total decay rate T{Hi,). Now the HQE is an operator product expansion 
(OPE) expressing r(i7(,) in terms of matrix elements of local AB = Q (B denotes the 
bottom number) operators, leading to an expansion of T{Hi,) in terms of AqcD/n^b- 
In the leading order of Aqcd /nT-b the decay rate of Hf, equals the decay rate of a free 
6-quark, unaffected by the light degrees of freedom of Hi,. Consequently, the lifetimes 
of all 6-flavored hadrons are the same at this order The dominant source of lifetime dif- 
ferences are weak interaction effects between the 6-quark and the light valence quark. 
They are depicted in Fig. |l] for the case of the B'^-B^ lifetime difference. The rela- 
tive size of these weak non-spectator effects to the leading free-quark decay is of order 
167r^(AQ(7£i/mf,)'^ — 0(5—10%). The measurement of lifetime differences among 
different 6-flavored hadrons therefore tests the HQE formalism at the third order in the 
expansion parameter. 

The optical theorem relates the total decay rate T{Hi,) to the self-energy of H),: 

- -^{H,\T\Hb). (1) 
Here we have introduced the transition operator: 

r = ImiJd^xT[H{x)H{0)] (2) 

with the effective \AB\ = 1 Hamiltonian H describing the T4^-mediated decay of the 
6 quark. The HQE amounts to an OPE applied to T which effectively integrates out 
the hard loop momenta (corresponding to the momenta of the final state quarks). We 
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decompose the result as 

T ^ [To + T2 + %\ [1 + 0{\lml)\ 

% = r^ + T'^ + r,™,,. (3) 

Here 7^ denotes the portion of T which is suppressed by a factor of 1/ with respect 
to To describing the free quark decay. The contributions to Ta from the weak interaction 
with the valence quark read 

+ (d ^ s) 

(^2 2 IT/ 12 

-jd ^ ^ ''I "''I [F'^Q" + i^lQ" + G'^T" + G|TJJ] . (4) 
67r 

Here Gf is the Fermi constant, is the bottom mass and the T^j's are elements of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix. The superscript q of the coefficients F'^, 
Fg, G^, G| refers to the cq intermediate state. The leading contributions to T" and 
T'^ are obtained from the left and right diagram in Fig. respectively. They involve 
the local dimension-6, Ai? = operators 

Q% = 6(1- 75)9 g(i + 75)fe, 

= b^^(l-^^)T-qq^^^{l-^^)T'^h, 
T| = &(l-75)r"<7g(l + 75)T'^&, (5) 

where T° is the generator of color SU(3). The Wilson coefficients i^" . . . Gg contain 
the physics from scales above rrib and are computed in perturbation theory. The remain- 
der Tsing in (^) involves additional dimension-6 operators, which are SU{'d) f singlets 
and do not contribute to the lifetime splitting within the B^) and S^) iso- 
doublets. In order to predict the widths T{B^) and r(S+) one needs to compute 
the hadronic matrix elements of the operators in (||). After using the isospin relation 
(B^IQ''^"!^^) = the matrix elements will enter T{B°^) - \(B+) in 

isospin-breaking combinations, which are conventionally parametrized as [jsj^] 

{B+\{Q--Q'')\B+) = /|M|i?i, {B+m^-Q%)\B+) = /|M|i32, 
(B+KT" - T'')\B+) = flMle,, {B+\iTli ~ T|)|i?+) = flM^ (6) 

Here fs and Mb are decay constant and mass of the B meson, respectively. In the vac- 
uum saturation approximation (VS A) onehcis Bi — 1, B2 = l+0{as{mb), Aqcd/ fTT-b) 
and £1^2 = 0. Corrections to the VSA results are of order 1/Nc, where iVc = 3 is the 
number of colors. 



2 Lifetime differences at next-to-leading order 



3 



We now find from (M and (p 



r(i30)-r(i3+) = 



127r 



iKdr-F" + \VcdrF^ - b. (?) 



Here we have introduced the shorthand notation 



/ F^{z) \ 



B 



( B, \ 

B2 

V ^2 J 



for q ~ d,u, c. 



(8) 



Since the hard loops involve the charm quark, the coefficient depends on the ratio 
z = rn^/ml. The minimal way to include QCD effects is the leading logarithmic 
approximation, which includes corrections of order a" ln"(mf,/Mvi/), n = 0, 1, ... in 
F'^ in (Q). The corresponding leading order (LO) calculation of the width difference 
in involves the diagrams in Fig. [l| Yet LO results are too crude for a precise 

calculation of lifetime differences. The heavy-quark masses in cannot be defined in 
a proper way and one faces a large dependence on unphysical renormalization scales. 
Furthermore, results for B12 and £12 from lattice gauge theory cannot be matched to 
the continuum theory in a meaningful way at LO. Finally, as pointed out in [^, at LO 
the coefficients F, Fs in are anomalously small. They multiply the large matrix 
elements parametrized by i?i,2, while the larger coefficients G, Gs come with the 
small hadronic parameters £12, rendering the LO prediction highly unstable. To cure 
these problems one must include the next-to-leading-order (NLO) QCD corrections of 
order <+iln"(mb/Afiv). 

The first calculation of a lifetime difference beyond the LO was performed for the 
Bs~Bd lifetime difference [^, where ©(a^) corrections were calculated in the SU(3)f 
limit neglecting certain terms of order z. In this limit only a few penguin effects play 
a role. A complete NLO computation has been carried out for the lifetime difference 
between the two mass eigenstates of the B^ meson in Q. In particular the correct 
treatment of infrared effects, which appear at intermediate steps of the calculation, has 
been worked out in [^. The recent computation in [|l|] is conceptually similar to the 
one in [||], except that the considered transition is AB = rather than AB — 2 and the 
quark masses in the final state are different. The NLO calculation of T{B^) — T{B^) 
involves the diagrams of Fig. ^. In [Q| the NLO corrections to T{B^) — T{B^) have 
been calculated for the limiting case z = 0. The corrections to this limit are of order 
z In z or roughly 20%. The first NLO calculation with the complete z dependence was 
presented in ^ and subsequently confirmed in [^. 



2 Lifetime differences at next-to-leading order 

The analytic expressions for the Wilson coefficients F^j''^^ — • • ■ C^' ^j"* — ^J"* 

are cumbersome functions of z involving dilogarithms. They depend on the renormal- 
ization scheme chosen for the AB ~ operators in (|^) and also on the renormalization 
scale /ig = 0{mi,) at which these operators are defined. These dependences prop- 
erly cancel between F'^ and B in physical observables like When our results for 
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Figure 2: WA contributions in the next-to-leading order of QCD. The PI diagrams are 
obtained by interchanging u and d and reversing the fermion flow of the u and d lines. 
The first line shows the radiative corrections to AS = operators, which are necessary 
for the proper infrared factorization. Not displayed are the diagrams Eg, E4 and D3_g 
which are obtained from the corresponding unprimed diagrams by left-right reflection 
and the reverse of the fermion flow. 



G 



^ ,^ are combined with some non-perturbative compu- 

tation of Bi, . . . e2, one has to make sure that the numerical values of these hadronic 
parameters correspond to the same renormalization scheme. Our scheme is defined by 
the use of dimensional regularization with MS [|[] subtraction, an anticommuting 75 
and a choice of evanescent operators preserving Fierz invariance at the loop level [|[]. 
Choosing further /io = rrib the desired lifetime ratio can be compactly written as 



1 = r(B+) 



0.0325 



cb\ 



0.04 



nib 



4.8 GeV 



B 



200 MeV 



(1.0±0.2)Bi + (0.1 ±0.1)^2 - (18.4 ± 0.9) ei + (4.0 ± 0.2) €3 1 . (9) 

Here t{B^) = 1.653 ps has been used in the overall factor. 

The hadronic parameters have been computed in [ [To| ] with quenched lattice QCD 
using the same renormaUzation scheme as in the present paper They read 

(Si, 52,61,62) = (1.10 ± 0.20, 0.79 ±0.10, -0.02 ±0.02, 0.03 ±0.01). (10) 

Inserting \Vcb\ = 0.040 ± 0.0016 from a CLEO analysis of inclusive semileptonic B 
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decays 1 13 1, the world average /s = (200±30) MeV from lattice calculations [|T^ and 
rrih = 4.8 ±0.1 GeV for the one-loop bottom pole mass into our NLO prediction 



reads 



compared to 



r{B+) 



1.053 ±0.016 ±0.017 



= 1.041 ±0.040 ±0.013. 



(11) 



(12) 



LO 



Here the first error is due to the errors on the coefficients and the hadronic parameters 
([To|), and the second error is the overall normalization uncertainty due to mb, \Vcb\ 
and Jb in (^. The Wilson coefficients also depend on the renormalization scale /ii at 
which the AB = 1 operators entering the diagrams in Figs, [l] and ^ are defined. This 
dependence stems from the truncation of the perturbation series and diminishes order- 
by-order in tts . The dependence on /ii is the dominant uncertainty of the LO prediction 
of the lifetime ratio. In Fig. ^ the /ii -dependence of the LO and NLO predictions for 
t{B^) / t{B^) — 1 is shown. The substantial reduction of scale dependence at NLO 



leads to the improvement in the NLO vs. LO results in (|1 1|),(|12[). Note that the NLO 

not be drawn from the old LO result. The result in ( |ll| ) 
measurement from the B factories Qll], n2] : 



calculation has firmly established that t{B'^) > r(_B"), a conclusion which could 

is compatible with recent 



'(S3) 



1.082 ± 0.026 ±0.012 
1.091 ±0.023 ±0.014 



(BaBar) 
(BELLE) 



The calculated Wilson coefficients can also be used to predict the lifetime splitting 
within the iso-doublet (SjJ ^ bus, ~ bds) with NLO precision. The corresponding 
LO diagrams are shown in Fig. Note that the role of T" and T'' is interchanged 
compared to the meson case with T" describing the Pauli interference effect. The 
lifetime difference between Ab ~ bud and is expected to be small, as in the case 
of i?" and Bj^, because it mainly stems from the small U-spin breaking effects in the 
matrix elements appearing at order 1 / m^. 

For Sft's the weak decay of the valence s-quark could be relevant: the decays 
AftTT", — > Ai,e~i'e and S° A;,7r° are triggered by s ^ u transitions and could 
affect the total rates at the 0{1%) level Once the lifetime measurements reach 
this accuracy, one should correct for this effect. To this end we define 



r(Sb) 



r(Sb) - r(Sb ^ AbX) = 



l-B(Sb^AbX) _ 1 



for Sb = Sb,Sj , 



(13) 



where i?(Sb AtX) is the branching ratio of the above-mentioned decay modes. 
Thus r(Sb) is the contribution from b ^ c transitions to the total decay rate. In 
contrast to the B meson system, the matrix elements of the four operators in (||) are not 
independent at the considered order in Aqcd / n^b- Since the light degrees of freedom 
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Figure 3: Dependence of t{B^)/t{B^) — 1 on fii/mt, for the central values of the 
input parameters and /io ~ rrii,. The solid (short-dashed) line shows the NLO (LO) 
result. The long-dashed line shows the NLO result in the approximation of i.e. z is 
set to zero in the NLO corrections. 



areinaspin-0 state, the matrix elements (5b|2(3| + Q*|S;,) and (Sb|2T| + T'^|S;,) are 
power-suppressed compared to those in (p^ (see e.g. [|[|]]). This, however, is not true 
in all renormalization schemes, in the MS scheme used by us 2Q| + Q'^ and 2T| + 
receive short-distance corrections, because hard gluons can resolve the heavy 5-quark 
mass. A priori one can choose the renormalization of e.g. Qg independently from Q'', 
so that (Sb|2(5| + Q^lSfe) = 0{AQCD/'iTLb) can only hold in certain renormalization 
schemes. This is also the case, if the operators are defined in heavy quark effective 
theory (HQET) rather than in full QCD. After properly taking into account these short- 
distance corrections, one can express the desired lifetime ratio solely in terms of two 
hadronic parameters defined as 



(S°|(T«-T'^)(Mo)|Sg) = /|MBMH,ir(Mo). (14) 
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ws 



Figure 4: Weak scattering (WS) and PI diagrams for S;, baryons in the leading order of 
QCD. They contribute to r(S|J) and r(S^), respectively. CKM-suppressed contribu- 
tions are not shown. 



Then one finds 



-1 = r(sg) [r(s-)-r(sg)] 

2 



= 0.59 



Vr 



cb\ 



U.SGeV/ 



/£ 



0.04/ V4.8GeV/ V200MeVy 1.5ps 

(0.04 ± 0.01) Li - (1.00 ± 0.04) i2 , (15) 



with Li = Lf" (/io = TOfe)- For the baryon case there is no reason to expect the color- 
octet matrix element to be much smaller than the color-singlet ones, so that the term 
with L2 will dominate the result. The hadronic parameters Li 2 have been analyzed in 
an exploratory study of lattice HQET [ [T^ for Ah baryons. Up to SU(3)f corrections, 
which are irrelevant in view of the other uncertainties, Lf" and L^' are equal. 



3 Conclusions 

Twenty years ago the ITEP group has developed the Heavy Quark Expansion, which 
allows to study inclusive decay rates of heavy hadrons in a model-free, QCD-based 
framework [|[]. The HQE expresses these decay rates as a series in both AQCD/m-b 
and as(TOfc). With the advent of precision measurements of lifetimes of 6-flavored 
hadrons at the B factories and the Tevatron correspondingly precise theory predictions 
are desirable. This requires the calculation of higher-order terms in the HQE. The 
inclusion of the corrections presented in this talk is in particular mandatory for any 
meaningful use of hadronic matrix elements computed in lattice gauge theory. The 
calculated QCD corrections to the WA and PI diagrams in Figs. |l|,^ allow to study the 
lifetime splitting within the (S^, B^) and (S°, Sj^) iso-doublets with NLO accuracy. 
It is gratifying that these corrections have been independently calulated by two groups 
finding agreement in their analytic expressions for the Wilson coefficients [|l|0]. 

Current lattice calculations, which are still in a relatively early stage in this case, 
yield, when combined with our calculations, t{B+)/t{B°i) = 1.053 ± 0.016 ± 0.017 
[see (11)]. The effects of unquenching and I /nib corrections are not included in the 
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error estimate, but the unquenching effects can well be sizable. A substantial improve- 
ment of the NLO calculation is the large reduction of perturbative uncertainty reflected 
in the scale dependence stemming from the AB — 1 operators. This scale dependence 
had been found to be very large at leading order, preventing even an unambiguous 
prediction of the sign of r(_B+)/r(i?°) — 1 up to now [^. 

At present the experimentally measured Af, lifetime falls short of t{B^) by roughly 
20% [|l7|], which has raised concerns about the applicability of the HQE to baryons. 
Unfortunately this interesting topic cannot yet be addressed at the NLO level for two 
reasons: First, r(Af,) /t{B^) receives contributions from the yet uncalculated SU(3)f- 
singlet portion Tgmg of the transition operator in (^. Second, the hadronic matrix 
elements entering t{Ai,)/t{B^) involve penguin (also called 'eye') contractions of 
the operators in (^), which are difficult to compute. These penguin contractions are 
contributions to the matrix elements in which the light q and q quark fields of the 
operator are contracted with each other, not with the hadron's valence quarks. 
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